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i-^ I We derive nonperturbative flow equations within an effective constituent quark 

JL , model for two quark flavors. Heat-kernel methods are employed for a renor- 

^ • malization group improved eff'ective potential. We study the evolution of the 



effective potential with respect to an infrared cutoff scale k at vanishing tem- 
perature. At the first stage we omit corrections coming from the anomalous 
dimension. This investigation is extrapolated to finite temperature, where 
we find a second order phase transition in the chiral limit at Tc ~ 130 
MeV. Due to a smooth decoupling of massive modes, we can directly link 
the low-temperature four-dimensional theory to the three-dimensional high- 
temperature theory and can determine universal critical exponents. 
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1 Introduction 



The linear a-model has become one of the most used working tools for the inves- 
tigation of the chiral phase transition in QCD. It shares an adequate realism with 
good practicality in applications. Besides the study of equilibrium phenomena, the 
(T-model has also been widely applied to calculations of the disoriented chiral con- 
densate. In a series of ground-breaking papers Wetterich and his group |^, Q have 
applied renormalization group methods to calculate the parameters of the a-model 
for all resolution scales. The fundamental idea is to follow the dynamics of the sys- 
tem by integrating out quantum fluctuations in infinitesimal intervals from a high 
momentum scale to the far infrared [0, ^. In practice these exact renormalization 
group flow equations have to be truncated in a similar way to Schwinger- Dyson 
equations. They allow to treat the critical fluctuations of the long range a- and 
TT-fields near the second order phase transition, where mean field methods or sub 
summations of the effective mass type are insufficient. This approach also has to 
make approximations on the form of the effective action. In the following paper we 
follow the spirit of this renormalization group approach, parameterizing the shape of 
the effective potential. In fact we only allow quartic and quadratic couplings in the 
effective mesonic potential. This approximation together with a novel heat-kernel 
infrared cutoff prescription allows us to derive new and very transparent formulae 
for the evolution equations. Since our cutoff function is different from previous work 
it allows to see effects of such a variant of the cutoff functions used in Wetterichs 
group |I|, ^. The main aim of our contribution is to give analytical insight into the 
physics inherent in the method of evolution equations. 

For high resolution, i.e. short distance processes fundamental QCD with quarks 
and gluons is the most efficient theory. A typical scale associated with perturba- 
tive QCD is A > 1.5 GeV. RHIC and LHC physics for secondary particles with 
Px > A will be dominated by such processes. In nucleus nucleus collisions, however, 
scattering at smaller momentum scale will be non negligible. Since the strength 
of the QCD interactions increases with decreasing momentum transfer, character- 
istic qq bound states will form and influence the dynamics at larger distances. A 
typical resolution where these processes start to become important is A^^^ ?« 1.0 
GeV. Below this scale the vacuum changes and acquires a quark condensate and/or 
meson condensate. In a recent paper on deep inelastic scattering ||^ it has been 
shown that a photon with varying resolution Q^ is a good physical probe to see the 
transition from partons to constituent quarks experimentally. Around Q^ = 1 GeV^ 
the behavior of the structure function F2 changes qualitatively as a function of Q^ 
indicating that nature knows about the phenomenon of chiral symmetry restoration. 
We think that our heat-kernel cutoff restricting the virtualities of the intermediate 
states corresponds to the resolution scale of the photon. We use the experimental 
indication for the transition at 1 GeV as input to our calculation. In the sponta- 
neously broken phase of chiral symmetry the constituent quark mass will increase 
with decreasing resolution and remain finite down to A ?« ^qcd- Around this scale 
the confining gluon configurations make themselves felt via confining forces. In the 
interval Aqcd < ^ < ^xSB the dynamics is governed by constituent quarks inter- 



acting via pions and a-mesons. One sees nicely in our approach how the different 
quantum ffuctuations from the cr-mesons and constituent quarks become unimpor- 
tant when the infrared scale parameter becomes smaller than the respective masses 
of these states. These modes then decouple from the further evolution, leaving the 
zero mass pions alone in the evolution. 

After having solved the evolution equations at zero temperature for reasonable 
starting values of the coupling constants, we pursue the evolution at finite temper- 
ature. Here the relevant parameter is the ratio of the temperature over the infrared 
scale parameter. Decoupling now sets in when the ratio of masses plus the Mat- 
subara frequencies over the infrared scale becomes large. At high temperatures the 
summation over Matsubara frequencies is dominated by the lowest mode, thereby 
reducing the dynamics to the corresponding three-dimensional field theory, which 
is the purely bosonic 0(4)-model. We find for the critical index j3 = 0.40 which is 
in good agreement with the results from Monte Carlo calculations 0, e-expansions 
nl andof [|, |. 



The outline of the paper is as follows: In section 2 we derive the evolution 
equations assuming a fixed parameterization of the effective potential. In section 
3 we give the explicit expressions of the differential equations for the broken and 
symmetric phase both for zero and finite temperature. Section 4 is devoted to a 
presentation and discussion of the results obtained after numerical integration of 
the evolution equations. 



2 Effective action with infrared cutoff 



In this section we show the derivation of the flow equation for the SU{2) x SU{2) 
constituent quark model (CQM). 

The partition function for the SU{2) x SU{2) model at zero temperature is given 

by 

Z[J = 0] = fvqVqVaVTfexp{-fd^x{CF + CB)} ■ (1) 

We omit external sources J and investigate the chiral limit. The Euclidean space 
lagrangian in d = 4 dimensions looks like 

Cp = q{x) {'jEdE + g (cr + iTTT-f^)) q{x) , (2) 



The T = parameters of the linear cr-model are fixed at the ultraviolet scale A = 
1.2 GeV in a similar way as in ref. |^. At this scale the quarks are massless partons 
and the cr-field has no vacuum expectation value. The mass squared m^ = (0.550 
GeV)^ is positive reflecting a symmetric ground state. The minimum of the effective 



potential U{a) lies at the origin. The other couplings are chosen as Aq = 40.0 and 
g = Mq/ /t, = 3.23, where g is the Yukawa coupling of the constituent quarks to 
the mesons and its value corresponds to a constituent quark mass Mg = 300 MeV. 
We do not evolve the quark Yukawa coupling in this introductory work. The result 
of the evolution will be largely insensitive to the exact starting value of Xq. We 
could also have chosen Aq = corresponding to an initial Lagrangian of the Nambu- 
Jona-Lasinio type without explicit mesonic interactions and non-propagating mesons 



12. 13 



Integration over the fermions yields formally a non-local determinant, which can 
be defined by a heat-kernel representation: 

Z[J = 0] = fVaViidet {-fEdE + gM{x)) exp{- f d'^xCs} 

= f VaVTfexp{--Tr log DD+ - fd^xCBJ (4) 

= / VaVif exp{— / d'^xCs + 

2 J1/A2 r J 

where we use the abbreviations M(x) = a{x) + ir7r(x)75 and D = •jEdE + gM{x). 
We remark that the lower boundary of the integral over the proper time r reflects 
the ultraviolet scale A which is fixed during the whole calculation. In fact, we 
will choose our infrared cutoff function in such a way that no ultraviolet divergences 
arise. Details concerning the heat-kernel representation and definitions can be found 
in ref. |^. 

In the following we are not interested in wave function renormalizations. This 
means that we set the wave function renormalization constant Zk = 1. We rewrite 
the meson fields in vectorial form with = (a, tt) and MM~^ = cfP'. We omit all 
derivatives in the heat-kernel expression and get for the partition function 



/ 



Z[J = 0] = 
VaVTieMl r — [d''xtr{x\e-<-^'^+'''^')\x)- fd^xCB} ■ (5) 

2 J1/A2 T J J 



In principle, the fermion determinant is part of the functional integration over 
the meson fields weighted by the meson lagrangian Cb- In the following, however, we 
will limit the integration over fermion and boson fluctuations to the one-loop level 
with the additional condition that the virtualities of the propagators in the loops 
are restricted to the interval [fc^, A^], where k"^ is the infrared and A^ the ultraviolet 
cutoff, i.e. we do not consider the effect of modified fermions on the meson dynamics 
in agreement with the one-loop approximation. The total effective action governing 

the dynamics of the slowly varying meson fields with virtualities < /c^ becomes a 
sum of fermionic and bosonic terms. 



Using a plane wave basis for the diagonal part of the heat-kernel the effective 
action for the fermions T^ can be defined by 

Here the remaining trace goes over color-, flavor- and spin-space. We leave out 
gradients of in the exponential and restrict the background field to low virtu- 
alities. The equivalent one-loop integration over the meson fields can be done by 

2 — * — * 

expanding the mesonic potential Vq = ^(0^) + ^(0^)^ around the slowly varying 
field configurations and one finds for the effective action for the bosons 

Here the trace goes over the (4 x 4) fluctuation matrices in {a, 7f)-space. Note the 
opposite sign of the bosonic and fermionic actions. The total effective action in 
one-loop approximation is given by the sum of both actions: 

r(0) = r^(0) + r^(0). (8) 

In order to implement in this effective action the infrared cutoff on virtualities 
one introduces a universal A;-dependent function fk{x = rk"^) into the proper time 
integrand. Note, that the proper time gives the inverse of the virtualities of the 
system. We mentioned this already in the discussion of the ultraviolet cutoff. The 
function fk has to satisfy some general conditions: Since the action F^ with infrared 
cutoff should tend to the effective action F at /c = 0, one must require that 

fm = 1 . (9) 

To suppress modes with small virtualities the function fk must also satisfy 

/fc(x-^oo)^0 . (10) 

This condition regularizes the infrared region of the effective action. In addition one 
needs a further condition on the first derivative of the cutoff function fk{x) in order 
to enforce that the flow equation for F^ is ultraviolet finite: 

fl{x) = -x'g{x) (11) 

with g{x) being a regular function in the vicinity of the origin. 

One possible choice for the function fk{x) which fulfills all these requirements is 

h{x) = e-^il+x + lx"") . (12) 



The effective action with virtuality cutoff function /^(tA;^) has the following form: 

Ffc[0] = fd'xVk{4>) (13) 



with 

V.,(^) = -l^"^A(rt=) /!«,,,, {e-<''-^)_,-n,'.»".^.j . (14) 

Because of the form of /^(rA;^) derivatives of the effective potential with respect 
to k are now infrared as well as ultraviolet regularized, so a further ultraviolet cutoff 
A is no longer necessary. In the a-model the second derivatives of the potential Vq 
evaluated at 0^ are given by 



= \o{<P'-^l)5^, + 2\oMJ (15) 

and the trace can explicitly be evaluated for A^ = 4 components. We find for the 
trace 

tr^ve""^ = 3e-^^«('^^'-'^io) + e-"^°^''^^'-'^'o) . (16) 



The first term on the r.h.s. of eq. (|16D corresponds to the eigenvalues associated 
with the three pions while the other term describes the cr-meson. Since the fermionic 
trace is diagonal in color-, flavor- and spin-spaces we get 

tre-^9^'^' = AN.Nfe-^^'^' . (17) 

Using these results we are able to calculate the derivatives of Vfc(0) which we need 
for the renormalization group flow equations. 



3 Renormalization group flow equations 

We want to study the evolution of different quantities with respect to the virtuahty 
or momentum scale k. Since we are dealing with the possibility that the chiral 
symmetry is spontaneously broken we have to distinguish between two phases. For 
large k we start in the symmetric phase which is defined by a vanishing vacuum 
expectation value (VEV) {(pk = 0) of the potential. Thus by means of the derivative 
with respect to the fields 

V := ^ (18) 

we can define the mass ml and the coupling constant A^ for the symmetric regime 

by 



T 



^"(0^ = 0. = O) . (19) 



Taking the derivative with respect to the scale k we find the following coupled sets 
of fiow equations for the symmetric phase (0| = 0): 

kdmi _ dnoi .^„. 

kdX, _ dV"{0) 

2^ " ^^^ • ^^^^ 

In the spontaneously broken region the VEV is finite {(pk 7^ 0) and the mass 
parameter m| tends to negative values. We prefer to parameterize the evolution of 
the potential in this region in terms of Xk and the minimum of the potential (pk, 
which is defined by 

V'M = . (22) 

This equation enables us to find the evolution of the minimum 0^. By again taking 
the derivative with respect to the scale k we find the fiow equation for 0| in the 
broken phase (0| 7^ 0). Note that we first evaluate the derivatives on the right 
hand sides for constant couplings with respect to k and then improve this one-loop 
expression by substituting the running couplings A^, ml and the VEV 0^ instead of 
of Ao, rnl and the VEV 0^^. This procedure corresponds to renormalization group 
improvement. Thus we find the general fiow equations for the broken phase 



kd<f>l k dV'{<Pl) kdV'{<Pl) 



2 dk 2V"{(f)l) dk Xk dk 

kdXk JV"{<Pl) ^^ ^^,„ 



(23) 



2 8k ^^T ' '' ""'"-' • ^^^) 

3.1 Evolution for T = 

The heat-kernel representation of the effective potential yields the explicit evolu- 
tion equations at T = 0. In the symmetric phase at zero temperature we use the 



equations (p!6| , |20| , ^lD and get 



k drnj ^ 3Xkk' 1 AN^g^k^ 

2 dk (47r)2(l + ^2/p)2+ (4^)2 ' ^ ) 

kdXk 12X1 1 8iV,(?' 



2 dk (47r)2(i+^2/p)3 (4^) 



(26) 



Note the Yukawa coupling g is fixed. In general we follow the evolution equations 
starting at fco = 1-2 GeV with the initial values Afc^ = Aq and m\^ = rriQ. We proceed 
from large k to small k integrating out more and more infrared modes. Near k = k^ 
the evolution of A^ is dominated by the bosonic term proportional to A| whereas the 
evolution of m^ responds to the fermion loop. We reach chiral symmetry breaking 
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Figure 1: The evolution of A^ with respect to the scale k. 



at the scale k = A^sb ~ 0.8 GeV, where we switch to the equations for the broken 
phase characterized by 0| 7^ 0. 

The squared mass of the cr-meson is then given by 2Afc0| and the quark mass by 
g(j)k and the evolution equations have the following form: 



k< 



se 



2 dk 

kdXk 
2~dk 



2(47r)2 
(47r)2 



1 + 



1 



a + 2A,02/P) 

3 

:i+2Afc02/p)3 



2^2 



ANc k^g 



(47r 



'9' 



[l + gm/k^ 



,(27) 



:281 



The first terms on the r.h.s. of the above equations are related to the bosonic while 
the last parts are connected to the fermionic contributions. The factor one in the 
first part comes from the three pions which are massless Goldstone bosons while the 
second part of the bosonic term contains the contribution of the a-meson. 

The (T-meson decouples from the evolution when the term proportional to 2\k(f)\/ 
fc2 becomes large, i.e. when the a-meson mass over the infrared parameter k is big. 
(cf. |]1|). The functions in squared brackets can be called threshold functions, since 
they describe the decoupling of the pions and cr-mesons from the evolution. With the 
heat-kernel method and the choice of the function /fc(T) these threshold functions 
can be obtained analytically. This is an advantage of our method compared to 
the momentum cutoff used in refs. [|l], Q. We remark that in the heat-kernel for 
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Figure 2: The evolution of the minimum of the effective potentiaL 



the effective potential the inverse fermion propagator enters quadratically in the 
combination DD~^ therefore we can use the same infrared cutoff function fk as in 
the bosonic integral without breaking chiral symmetry. If we want to evaluate the 
running coupling more care has to be used to regulate the fermion integration. One 
also recognizes the signs of the bosonic and fermionic contributions to the /3-function. 
The bosons lead to an infrared stable (ultraviolet unstable) coupling, whereas the 
fermions counteract this tendency. Going from high k to low k one sees that the 
mesonic self-interaction A^ will balance at intermediate values of k, whereas in the 
far infrared the boson term wins (cf. figure]^). At k = k^sB ~ 800 MeV eqs.(p6[) and 
(^) become identical, also visible in figure |l] and the /3-function of Xk is continuous. 
The vacuum expectation value (pk stabilizes at small values of k and the evolution 
ends with linifc^o 4'k = fn ■ When the heavy particles, the sigma and quarks, have 
decoupled, the change of the vacuum expectation value becomes proportional to k 
which vanishes (see figure ||). The infrared k^sB scale found in this calculation is 
somewhat smaller than the resolution Q of the photon found in an analysis of deep 
inelastic scattering 0] for the transition from the quark as a parton to the massive 
constituent quark. One has to wait for a more sophisticated calculation with running 
Yukawa coupling g^ and wave function renormalization Zk [^ which may improve 
the already astonishing agreement between phenomenology and the field theoretic 
model. 



The above equations are the evolution equations for the running of the expecta- 
tion value (pl and quartic coupling Xk with the infrared scale k at zero temperature. 



3.2 Evolution for finite T 

At finite temperature we integrate the momenta in equation (|T^) , by splitting the 
zero component from the three-dimensional spatial components and convert the 
integration over go into a summation over Matsubara frequencies u"^ = An'^n^T'^ for 
the bosons and u^ = {2n + l)^7r^T^ for the quarks. 

In the symmetric phase the corresponding equations are written in terms of the 
positive mass parameter m\ and A^: 



k drni 



2 dk 



k dxl 

2 dk 
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Figure 3: The bosonic threshold functions for different mass pa- 
rameters y = 0,2,4 as function of T/k. The dashed line is the 
function ISvr/S ■ T/k which demonstrate the linear behavior of the 
threshold function for large T/k. 
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In the broken phase one finds: 



2 dk 



3k^ 



2(47r)2 

_ AN, g^e 



SvrT ~ 
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2 A:„-^|(i+^2/p)./^ ■ (1 + (^2 + 2An0D^) A^) 
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(31) 



il + ^l/k'Y^' (i + K + 2An0D^)A^) 
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8N, 



(47r 
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ISvrT ^ 



;i + (^^ + ^^(0iP)A^) 



\7/2 



(32) 
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Figure 4: The fermionic threshold functions for different mass 
parameters y = 0, 2, 4 as function of T/k. 



One sees that the form of the equations is unchanged concerning the dependence 
on the couphng constants. Only the threshold functions given in square brackets are 
different at finite temperature. Due to the three dimensional momentum integrations 
fractional powers arise in the threshold functions. 

In the limit of low temperatures T/k ^ we regain our old expressions at T = 0, 
i.e. we find the relations: 
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SvrT ^ 1 1 



and also for the other threshold functions 

These equations guarantee the right matching of the finite temperature equations 
to the zero temperature equations, i.e. in the limit T — > the set of equations for 
finite and zero temperature become identical. For large ratios T /k the bosonic 
threshold functions, given in the square brackets of eqs. (^)-(^), increase linearly 
in T /k. (cf. figure |]). 

This increase is due to the n = Matsubara mode in the frequency sum. In the 
large temperature limit the (3 + l)-dimensional system reduces to a 3-dimensional 
system. We plot in figure |^ the threshold functions for the bosons and in figure ^ 
the threshold functions for the fermions. The dimensional reduction can also be 
seen in in the critical behavior, as will be shown in the next section. The fermionic 
threshold functions decrease with T /k and run to zero. There also exists a stable 
plateau in the vicinity of the origin for the Fermi-case as in the ref. . 



4 Numerical results and discussion 



The procedure to obtain finite temperature results is now the following. For each 
fixed temperature T we solve the evolution equations as functions of k with the same 
starting parameters as the T = theory. The underlying idea is that at the large 
ultraviolet scale the finite temperature does not modify the effective theory, since 
the finite temperature only modifies the boundary in the imaginary time direction. 
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Table 1: Different initial values compared with the corresponding 
critical temperatures. 



At the initial cutoff scale the momenta at the ultraviolet cutoff are anyhow so 
large that the theory is not affected by temperatures 27rT < 1.2 GeV. Of course, 

12 



for higher temperatures also the input would have to be modified. With increasing 
temperature the mass parameter ruk decreases slightly more slowly towards the 
condensation point (cf. figure ^, but the main effect of the finite temperature occurs 
below k ~ 800 MeV. 
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Figure 5: VEV (pk and mass ruk as function of k for different 
temperatures. (Upper line T = 11 MeV, middle line T = 111 MeV 
and bottom line T = 131 MeV). 



The threshold functions contain an additional damping due to the Matsubara 
frequency. Below the condensation point the boson condensate (pl reaches less high 
values at finite temperatures than before at T = 0. For all temperatures it decreases 
again with fc — > 0, which is due to the pion loop. After optimizing the Runge-Kutta 
code a critical temperature 



Tc ^ 131, 5 MeV 



(35) 



was found. 



This temperature would be compatible with that of the Wetterich group if /^r = 
93 MeV had been used in [^ even in the chiral limit [^. They obtained a critical 
temperature T^ ^ 100 MeV by using non-analytic threshold functions and wave 
function and coupling constant renormalization. Both approaches share the linear 
sigma model with free quarks even at low temperature. In table |l| other possible 
initial values for Xkg and rrikQ at the ultraviolet scale are shown which would give the 
same f.^ value for vanishing temperature. The corresponding critical temperature 
and k^sB scale are almost constant and are correlated like 



27rTc^ 



"XSB 



(36) 
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Figure 6: The order parameter as function of T. 



In figure |]we show the normahzed order parameter 0fc=o(7')/0fc=o(O) as a func- 
tion of temperature. It follows the result of chiral perturbation theory until T ^ 35 
MeV, then it deviates because of the stronger effects of the quark loops. In chiral 
perturbation theory the temperature dependence of the light quark condensate for 
massless quarks, which is plotted in figure 0, is given by the expression 



(gg)c 



1- 



T' 



rpA 



8/2 384/4 



-t In ^ + OiT^ 

288/6 T ^ 



(37) 



with Ag = 470±110 MeV |TH]. The purely mesonic description of the phase transition 



with a finite number of mesons becomes inadequate since more and more mesons 
become important. In the vicinity of Tc we also extrapolated in figure ^ the scaling 
behavior of the order parameter with our critical exponent 13 ~ 0.40. Still in a model 
without confinement of quarks the fermion loop may be overestimated cf. [0. The 
renormalization group allows to include the long range correlations properly near 
the critical point. In fact it is possible to see scale invariance by analyzing the power 
law behavior of the order parameter near T^. We plot log{(f)k) versus log{Tf. — T) in 
figure 0. The data points lie on a linear curve given by 



Q.mog{T^ - T) + (1 - 0.40)/o^(Tc) - 0.30 . 



(3^ 



This yields a critical exponent (5 ~ 0.40 which is in good agreement with lattice 
results for the 0(4)-theory in three dimensions [§]. So in spite of the finite temper- 
ature the limit T /k —>■ cxd enforces the dimensional reduction. One knows, that the 
linear sigma model lies in the same universality class as the 0(4)-theory, since the 
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fermions are not contributing to the critical fluctuations near Tc IQ- To compare 
the result with the mean field value (3 = 0.5 we also plotted this slope in the figure |^. 
Further investigations must study how far away from the critical point the renor- 
malization group improvement is important for the behavior of the order parameter. 
An especially interesting point is to find the window for the critical dynamics. In the 
superconducting phase transition only the mean field behavior is experimentally rel- 
evant. From the figure ^it seems that in chiral QCD the critical behavior infiuences 
a wider range of temperatures. In real QCD, however, the finite quark masses spoil 
the second order phase transition. Also gluon effects will change the behavior of the 
pressure in comparison with the effective linear a-model. Numerical simulations of 
lattice QCD still show an inconclusive critical behavior of the order parameter. In 
lattice QCD the main problem is the treatment of light particles leading to large 
correlation lengths which are larger than the lattice size. 

The presented method of heat-kernel regularization opens the way for more re- 
fined calculations also including coupling constant and wave function renormaliza- 
tion. Here the more delicate problem of introducing a cutoff which does not violate 
chiral symmetry has to be dealt with. Besides comparing the T = theory to 
the behavior of deep inelastic scattering we plan to investigate the spectrum of the 
Dirac operator [0. Since our method is a multi-scale analysis of the effective theory 
we think that the Dirac eigenmodes up to fc = 1.2 GeV can be checked and com- 
pared to lattice calculations. Random matrix theory has been extremely successful 
to predict the smallest eigenvalues and correlations of nearest neighbor levels. Our 
method should give the broad behavior of the Dirac spectrum at T = and for finite 
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temperatures. In our method the real dynamics of QCD at low resolution seems to 
be captured quite well. So we are optimistic about further prospects. These include 
the discussion of three flavor dynamics and finite quark mass effects [|18]. Due to 



the simplicity of our cutoff function we expect also a considerable simplification 
in the set up and solution of three flavor dynamics. For nuclear physicist remains 
the ever challenging field of finite baryon density where analytical methods like our 
renormalization group equations are very promising to model the transition of the 
nuclear many body system to the quark many body system. 
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